Ab-initio thermodynamics of deposition growth: surface terminations of CVD 

titanium carbide and nitride 
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We present a calculational method to predict terminations of growing or as-deposited surfaces 
as a function of the deposition conditions. Such characterizations are valueable for understanding 
catalysis and growth phenomena. The method combines ab-initio density functional theory (DFT) 
calculations and experimental thermodynamical data with a rate-equations description of partial 
pressures in the reaction chamber. The use of rate equations enables a complete description of a 
complex gas environment in terms of a few, (experimentally accessible) parameters. The predictions 
are based on comparisons between free energies of reaction associated with the formation of surfaces 
with different terminations. The method has an intrinsic non-equilibrium character. In the limit 
of dynamic equilibrium (with equal chemical potential in the surface and the gas phase) we find 
that the predictions of the method coincide with those of standard ab-initio based equilibrium 
thermodynamics. We illustrate the method for chemical vapor deposition (CVD) of TiC(lll) and 
TiN(lll), and find that the emerging termination can be controlled both by the environment and 
the growth rate. 

PACS numbers: 68.55.A-, 81.15.-z, 68.35.Md, 05.70.Ln 
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I. INTRODUCTION 

The structure and chemical composition of surfaces 
play a fundamental role in many industrial applications. 
Identifying and ultimately learning to control surface ter- 
minations is a key issue in modern materials design. 
In heterogenous catalysisA chemical reactions typically 
take place on the solid surface of the catalyst. Surfaces 
with different chemical compositions, for example, dif- 
ferent surface terminations, support different chemical 
reactions.— - — This is true even if the exposed surfaces 
possess identical crystallographic orientations. Similarly, 
surface terminations strongly influence nucleation and 
growth processes. Different surface terminations do, in 
general, favor adsorption of different atomic species^ - — 
and therefore affect subsequent growth of multicompo- 
nent materials. 

Atomistic modeling^ - — is highly valuable for design of 
functional materials and surfaces. It serves as an impor- 
tant complement to and extension of experimental char- 
acterizations. It also provides predictive power and thus 
an opportunity to accelerate innovation^ Coupling DFT 
to either thermodynamic" - — or kinetic modelings- 
has proven extremely useful for interpreting and com- 
plementing experimental techniques in characterization 
of surfaces^ thin films, 24 and interfaces^ Such model- 
ing is often well-suited for descriptions of structures that 
are fabricated or investigated under well-controlled con- 
ditions, such as ultra-high vacuum (UHV) or molecular 
beam epitaxy^ (MBE). 

Industrial mass production of materials often employ 
more efficient but less-controlled deposition methods, 
for example, chemical vapor deposition (CVD).— Such 
methods give rise to compositions and structures that 
need not be stable in a thermodynamical sense under am- 



bient or ultra-high vacuum (UHV) conditions. The fast 
deposition causes a lot more complexity in the growth 
proccess than is found in MBE. For example, it causes 
cluster formation and formation of intermediates in the 
gas phase or on the substrate and these structures de- 
fine barriers for subsequent deposition events. Relevant 
experimental characterizations of surface terminations in 
such cases are difficult and require in-situ measurements. 
Atomistic, ab initio based modeling is consequently of 
even higher value for as-deposited structures. However, 
a simple scale-up of the kinetic modeling 19-22 faces enor- 
mous practical problems because complex, unknown, and 
constantly evolving transition states determine the kinet- 
ics of the deposition. 

In this paper, we propose an ab-initio thermodynam- 
ics method for deposition growth (hereafter referred to as 
AIT-DG). Our method combines DFT calculations with 
thermodynamic concepts and rate equations. The latter 
effectively describe supply and exhaust of gases to and 
from a reaction chamber as well as deposition processes 
inside. Steady-state solutions of these rate equations en- 
able us to determine Gibbs free energies of reaction 2 ^ of 
deposition processes that lead to different surface termi- 
nations. The free energies of reaction allow us to predict 
surface terminations as a function of the growth environ- 
ment. 

We emphasize that the AIT-DG method is fundamen- 
tally different from the ab-initio thermodynamics method 
proposed in Refs. [HI and [H also adapted to interfaces^ 
and used by many others. Such methods are designed 
to predict terminations of surfaces that are in equilib- 
rium with a given gas environment (and are hereafter re- 
ferred to as AIT-SE, ab-initio thermodynamics with sur- 
face equilibrium). In contrast, our AIT-DG formulation 
possesses a non-equilibrium character, but formally con- 



tains the AIT-SE method in the regime, sometimes called 
dynamic equilibrium^ 2 - where one can assume equality of 
the chemical potentials in the surface and in the gas. 

We illustrate the AIT-DG method for the CVD 
TiX(lll) (X = C or N). These belong to the simplest 
class of materials where different surface terminations can 
arise. At the same time there is a high industrial interest 
in these materials. TiN is widely used as diffusion bar- 
rier for Al-based interconnects^ and finds application as 
ohmic contact in GaN semiconductor technology. 31 TiC 
can be used as substrate for growth of other carbidic ma- 
terials such as SiC or graphene^ 2 - and its (111) surface 
also shows a potential as catalyst^ Moreover, TiC and 
TiN both are commonly used as TiX/alumina multilayer 
coatings in cutting tool industry^ 

The paper is organized as follows. SectionlPIlfocuses on 
the discussion of surface terminations of binary materials. 
In Section Hm we present our AIT-DG method in detail, 
exemplified for CVD of TiX, and also considering the 
equilibrium limit of our method. Section IIVI summarizes 
the computational method that we use for all ab-initio 
calculations. We present our results in Section [V] and 
discuss these as well as the method itself in Section I VII 
Section [VIII contains a summary and our conclusions. 
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FIG. 1. Different surface terminations for a film deposited on 
a substrate. The schematics shows the case for a binary ma- 
terial AB formed of alternating A and B layers, represented 
by differently colored balls. While the interfacial composition 
largely is determined by the substrate properties, the compo- 
sition of the surface layer strongly depends on the deposition 
environment, possibly leading to an AB film that slightly de- 
viates from full stoichiometry. 



A more involved analysis is necessary to understand what 
termination will emerge in deposition growth. 



II. SURFACE TERMINATIONS 

Figure Q] shows a schematics of thin films of a binary 
material AB (deposited on a substrate), with the pan- 
els AB:A and AB:B differing essentially only by their 
surface terminations. Here and in the following, AB:A 
(AB:B) denotes an AB surface with A-termination (B- 
tcrmination) . The figure illustrates the character of our 
modeling where the difference in surface terminations is 
represented by films of different thicknesses; formally, the 
left panel represents a stoichiometric film grown to 2n 
layers, while the right panel shows a nonstoichiometric 
film grown to 2n + 1 layers^ This case represents the 
simplest case of surfaces that can have more than one 
termination, and the question arises which one will be 
created in a specific fabrication process. 

Well-defined surface terminations can be created in 
a variety of different deposition methods, including the 
well-controlled MBE and less-controlled methods, such as 
CVD. MBE employs UHV and allows, in principle, one 
atom (or dimer) to be deposited at a time. An abrupt ter- 
mination of the growth process therefore provides good 
explicit control over surface coverage and termination. 

Less-controlled deposition methods typically yield 
much higher growth rates and are therefore more pop- 
ular for large-scale manufacturing. In CVD, a substrate 
is exposed to an inflowing supply-gas mixture at high 
temperatures and at, compared to MBE, high pressures. 
Since there is no explicit control on the number of indi- 
vidual atoms or molecules in the gas environment or on 
the surface, one cannot simply end the deposition pro- 
cess so as to create a pre-specified surface termination. 



II. 1. Materials example: TiC & TiN 

TiX(lll) (X = C or N) are important, specific, exam- 
ples where surface terminations have industrial relevance. 
Specifically, alternating deposition of TiX and AI2O3 by 
CVD 36,37 is commonly used to produce wear-resistant 
coatings on cemented-carbide cutting tools.— 

The top panel of figure [3] illustrates the CVD process 
used for TiX growth. A TiCl 4 -CH 4 -H 2 -HCl (for TiCSi) 
or TiCl 4 -N 2 -H 2 (for TiN 37 ) supply gas mixture is injected 
into a CVD chamber via valve 1. The chamber is kept at 
a high constant temperature and at a constant pressure. 
The overall reactions that lead to the deposition of TiX 
are 

TiCl 4 + CH 4 TiC + 4HC1, (la) 
TiCl 4 + ^N 2 + 2H 2 -> TiN + 4HC1. (lb) 

Unused supply gas as well as reaction products are ex- 
hausted from the chamber via valve 2. 

In TiX/alumina multilayers, the TiX(lll) surface ter- 
minations play a crucial role for the nucleation of alu- 
mina and consequently for the nature of the TiX/alumina 
interface in the coatings. In particular, Ti-terminated 
surfaces indicate an interface composition of the type 
TiX:Ti/0:Al203, whereas X-terminated surfaces favor 
a TiX:X/Al:Al 2 3 interface^^ Here, (0,A1):A1 2 3 
specifies the dominant chemical species (O or Al) of alu- 
mina at the interface. 

To the best of our knowledge, experimental studies on 
the TiX(lll) surface termination exist only for TiC^Sr— 
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FIG. 2. Schematics of chemical vapor deposition (CVD) of 
TiX and of the kinetics in the formation of excess layers. The 
upper panel gives a principle sketch of a CVD apparatus. A 
supply gas mixture that contains several different gases with 
different concentrations Ci is injected into the CVD chamber 
via valve 1 at a rate i? aU p. The chamber is kept at a constant 
temperature T and pressure p. Inside the chamber the gases 
react and form TiX. Reaction products and unused supply 
gas are exhausted at a rate i? ox h via valve 2. The set of 
lower panels illustrates kinetic effects that specify the surface 
termination. A Ti (X) surface layer can be deposited from a 
Ti (X) carrying gas, here shown for TiCU. As indicated, the 
process involves a set of unknown complex transition states 
(TS). 



but not for TiN. They all report a Ti-terminated surface. 
However, in all these experiments, the TiC(lll) surface 
was prepared by high-temperature annealing (and associ- 
ated selective evaporation) under UHV conditions. The 
observed Ti-termination of the TiC(lll) surface under 
such experimental conditions can be understood from a 
theoretical argument due to Tan et ai.^ estimating the 
selective high-temperature evaporation from the total- 
energy cost of pulling atoms out of material. 

A predictive, atomistic theory of the surface termina- 
tion of as-deposited TiX(lll) as a function of the growth 
environment remains of central value. The characteriza- 
tions and predictions are important, for example, to de- 



fine a proper starting point for understanding atomistic 
processes in the fabrication of TiX/alumina multilayer 
coatings^ 



II. 2. Kinetics of deposition processes 

The bottom panel of Fig. [5] shows a schematics of the 
kinetics leading to a surface system with an excess Ti 
layer. The initial state is represented by a TiCU molecule 
above the X-terminated TiX(lll) surface. In the final 
state, a Ti atom (or a layer after a number of events) is 
deposited. Additional molecules that enter or leave the 
reaction and produce the final state are not shown. We 
also sketch the energy landscape as a function of a generic 
reaction coordinate. We note that only the relative en- 
ergy positions of the initial and final states are shown in 
correspondence with the actual results of our modeling; 
the nature and energy of the relevant transition states 
are unknown. 

Surface morphologies obtained in MBE growth can of- 
ten by understood from kinetic Monte Carlo (kMC) sim- 
ulations, using DFT calculations to determine process 
barriers*^— Such simulations follow the surface mor- 
phology over some time by allowing for a number of 
events to take place at random. Conceptually it is simple 
to generalize the kMC approach to CVD, but in practice 
such a generalization will be extremely challenging. The 
chemical species supplied in CVD are usually more com- 
plex than single atoms or dimers used in MBE. Worse, 
the reactants will not, in general, simply dissociate on the 
surface but enter into complicated reactions with other 
species and form clusters. These reactions are not limited 
to the surface and may lead to a large variety of inter- 
mediates. While adaptive intelligent kMC exists and is 
being developed^ the complexity in CVD is enormous. 

At the same time, the high deposition temperature in 
CVD motivates a shift away from explicit kinetic mod- 
eling. We therefore develop a thermodynamic descrip- 
tion that focuses entirely on initial and final states, re- 
taining the simplicity of the analysis for selective evapo- 
ration in IJHVJ^ A nonequilibrium thermodynamic de- 
scription becomes applicable when high temperature fa- 
cilitates structural reorganization and permits the grow- 
ing system to reach and sample the set of possible final 
states. 



II. 3. Excess-layer deposition 

Figure [3] illustrates possible processes that lead to de- 
position of single Ti and C layers for TiC. We determine 
a preference for the surface termination by comparison 
of the Gibbs free energy of reaction for such excess-layer 
structures with those of the stoichiometric TiX. The re- 
actions listed in (JTJ) only describe deposition that leads to 
TiX with the full stoichiometry of the bulk. We supple- 
ment that description by a thermodynamic characteriza- 



tion of dominant reaction pathways leading to formation 
of excess layers for both TiC and TiN systems. 

For excess-layer reactions, we consider the following 
reactions for Ti deposition on TiC and TiN, 



TiC:C - 
TiN:N ■ 



TiCl 4 
■ TiCl 4 



2H 2 
2H 2 



TiC:C/Ti- 
TiN:N/Ti ■ 



4HC1, 
4HC1. 



(2a) 
(2b) 



Here, TiC:C/Ti (TiN:N/Ti) identifies a system with 
an additional Ti layer on the C-terminated TiC (N- 
terminated TiN) surface. Additional X layers may be 
deposited as 

TiC:Ti + CH 4 -> TiC:Ti/C + 2H 2 , (3a) 

TiN:Ti + -N 2 -> TiN:Ti/N. (3b) 

Since deposition of Ti is relevant only on an X-terminated 
surface and deposition of X only on the Ti-terminated 
surface we suppress the specification of the original sur- 
face termination in the following. For example, TiC/Ti 
will be used instead of TiC:C/Ti. 



III. AB-INITIO THERMODYNAMICS OF 
DEPOSITION GROWTH 

Our thermodynamic method to predict terminations 
of growing surfaces describes phenomena that lie between 
the static equilibrium limit and the non-equilibrium limit 
that requires a fully kinetic theory. We assume an exper- 
imental setup similar to that illustrated in the top panel 
of figure [21 



III.l. Computational strategy and assumptions 



We use the difference in Gibbs free energies of reaction 



AG, 



AB 



(4) 



as a predictor for the prevalence of growing either an A- 
terminated (AG^ B < 0) or a B-terminated (AG^ B > 0) 

surface of a binary material AB. Here, G^ B ^ is the free 
energy of reaction associated with the formation of an ex- 
cess A (B) layer. Our choice of predictor follows naturally 
from an analogy with chemical reaction theor y 28 ' 29 and 
from considering surface transformations. The following 
coupled set of reactions change an A- (B-) terminated 
AB surface into a B- (A-) terminated surface, 



I: B 



React/ f± A + Prod/ 



II : A + React// ^ B + Prod//. 



(5a) 
(5b) 
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FIG. 3. Formation of excess layers and associated free ener- 
gies of reaction G T . The top panel illustrates the formation 
of a C excess layer from CH4 on a Ti-terminated surface. 
The bottom panel shows the formation of a Ti excess layer 
from TiCU and H2 on a C-terminated surface. From a kinetic 
perspective, formation of specific clusters and presence of H2 
may catalyze the process shown in the upper panel; thermo- 
dynamically, such catalytic processes do not, however, affect 
the free energy of reaction. 



M f b\ depend on the deposition environment (reactants 



{/,*>} 

React{/ //} and products Prod{/ m) 
solution Pa/Pb = (T' f 
proximeted by 



. The steady-state 
Ty 7 ) is closely ap- 



Pa/Pb « exp (-/3AG^ B /2 



(0) 



where /3 denotes the inverse temperature. This evalua- 
tion ([6]) becomes exact in the limit of dynamic equilib- 
rium (identified by G^ + G B = 0) where AG^ B reduces 
to a difference in surface free energies. Further details 
on the nature of the strategy and the applicability to 
noncquilibrium conditions are given in the appendix. 

The assumptions that facilitate the evaluation of free 
energies of reaction are: 



The probability for growing either an A- or B-terminated 
surface is described by rate equations where the rates 



1. The system composed of the growing surface and 
the gas environment settles into a steady state, 
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specified by the gas supply, deposition, and gas ex- 
haust rates. We do not a priory assume that the 
gases and the surface are in mutual equilibrium, 
nor do we generally expect such an equilibrium to 
emerge. 

2. Each of the gases that take part in reactions is fully 
thermalized at the steady-state partial pressure so 
that we can associate a chemical potential to them. 

3. The overall pressure is low enough so that a de- 
scription of all gases in terms of the ideal-gas ap- 
proximation applies. 

4. The temperature of the system is high enough so 
that we can neglect kinetic barriers. 

The computational strategy and the first assumption of 
a steady state is a novel feature of our method. We derive 
the steady state from a rate equation that describes the 
supply and exhaust of gases into and from the reaction 
chamber as well as the deposition inside the chamber. 

The other three assumptions are identical to those of 
the AIT-SE method. The AIT-SE describes a system 
in a static equilibrium at high temperatures. The AIT- 
SE might also describe systems in dynamic equilibrium. 
Dynamic equilibrium^ here means that the free energy 
of reaction (derived and calculated below) exactly equals 
zero. This condition also applies in static equilibrium but 
in dynamic equilibrium the zero- value of the free energy 
of reaction arises from the steady-state concentrations of 
reactants and reaction products. A (virtual) drop in the 
reaction rate readily leads to a drop of the concentrations 
of reaction products but not of the reactants (due to the 
supply and exhaust), so that the dynamics of the reaction 
is restored. 

We find that as the steady state approaches dynamic 
equilibrium, the predictions of AIT-DG coincide with 
those of the AIT-SE. However, even if dynamic equi- 
librium is reached, our description in terms of a rate 
equation will be necessary in order to determine chemical 
potentials of the gas-phase constituents from their indi- 
vidual pressures during the process. We also emphasize 
that, during growth, the steady state does not necessarily 
settle into a state of dynamic equilibrium. 



along the [111] direction. In particular, it possesses one 
Ti-terminated surface and one X-terminated surface. An 
excess Ti layer is deposited on the X-terminated side of 
the slab, while an excess X layer is be deposited on the 
Ti-terminated side. 

We introduce the chemical potentials /tt^f, (J-Q g , and 
/i^ Ig for Ti, C, and N atoms in the gases. The set of /i^ lg 
can be conveniently related to differences between several 
molecular chemical potentials (see below). The chemical 
potentials of Ti, C, and N in the gases do not in general 
equal the respective chemical potentials /iTi, (J-C, and 
in the solid bulk. For the latter the free energy per bulk 
unit^ 5Tix, always specifies material stability, 

MTi + Mx — 9Tix- (7) 
Growth of TiX, on the other hand, is characterized by 

5TiX</4 S + A4 iS - (8) 

The free energy of reaction^ for the deposition of one 
unit of bulk is defined as 

TiX _ aig aig , Q \ 

We also introduce the free energy of reaction G r which 
corresponds to the gain in energy per surface unit in the 
deposition of one excess layer. For deposition of excess 
Ti, G r can then be written as^ 

GjiC/Ti = Gt . c/t . _ Gt . c _ ^aig (1Qa) 

G™/ Ti = G TiN/Ti - G TiN - 4f- (10b) 

Here, Gxix is the free energy of a stoichiometric, non- 
symmetric TiX slab and Gxix/Ti denotes the free energy 
of the system that consists of the same stoichiometric 
TiX slab and an additional Ti excess layer (adsorbed on 
the X-terminated side of the slab) . For excess X, the free 
energy of reaction is given by 

G™/ c = G TiC/c -G TiC -Mc is (Ha) 
G™/ N = G TiN/N - G TiN - MN iS - (Hb) 



III. 2. Excess layers and their free energies of 
reaction 

Figure [3J illustrates the formation of Ti and C excess 
layers on TiC(lll), corresponding to the reactions in (|2"aj) 
and (|3a|) . Similarly, we assume that Ti and N excess 
layers can be deposited on TiN(lll) according to the 
reactions in (f2b|) and (f3b|) . 

We represent the TiX(lll) surface before deposition 
of additional layers by a stoichiometric slab. Upon de- 
position of an additional layer the slab becomes non- 
stoichiometric, which motivates the terminology of an ex- 
cess layer. A stoichiometric TiX(lll) slab is asymmetric 



With these definitions, a negative value of G r implies 
that the deposition of the excess layer is thermodynami- 
cally favorable and that the reaction proceeds in the di- 
rection of the arrow in one of the listings and ([3]). 
A positive value of G r implies that the deposition of the 
corresponding excess layer is thermodynamically unfa- 
vorable and that the reaction proceeds in the opposite 
direction. 

The atom-in-gas chemical potentials fi^ ls depend on 
the composition of the environment and can therefore 
be controled by the gas flow. For the case of CVD of 
TiX with the here-considered supply gas and the reac- 
tions listed in @ and ([3]) , we assign the atomic chemical 
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potentials in the gas as 

^Xi = MTiCU + 2/iH 2 — 4^HC1 



III. 4. Rate equation and steady-state solution 



aig 



aig 



MCH 4 _ 2/in 2 
1 



(12a) 
(12b) 

(12c) 



We note that, in general, the deposition of excess Ti 
followed by the deposition of excess X (or vice versa) 
equals the free energy of reaction for the deposition of 
one unit of bulk, 



Gj iX / Ti + Gj ix / X = ffr TiX 



(13) 



If the solid-gas system approaches equilibrium (mediated 
through the surface), we have gJ lX — 0, so that the re- 
action energy gain is simply 



G TiX/Ti + G TiX/X = Q _ 



(14) 



We emphasize that the definitions (|III.2|) and (jTTJ) do 
not make any assumption about equilibrium. On the 
contrary, by calculating (1131) as a function of deposition 
pressure and temperature, we can predict whether the 
process is in or out of equilibrium, and in which direction 
it proceeds. 



III. 3. Gibbs free energies and chemical potentials. 

Calculating free energies of reactions requires simple 
approximations to evaluate free energies of surface sys- 
tems and chemical potentials of gas-phase species. We 
base these approximations on experimentally available 
thcrmochemical data. 

For the gases, we employ the ideal gas approximation. 
This approximation in justified for low pressures which 
is typically the case in CVD of TiX. It also allows us 
to express the fit in terms of temperature T and partial 
pressure pi , 

tii(T, Pi ) = et + A°(T) + k B Tln( Pt /p°). (15) 

Here ej is the DFT total energy of the gas phase species 
(molecule), and k B is the Boltzmann constant. A?(T) is 
the temperature dependence of fa at a fixed pressure p , 
and related to enthalpy and entropy differences^ The 
latter are available for p = 1 atm for many molecules in 
thermochemical tables such as Ref. l45l 

For surface systems, we approximate the Gibbs free 
energy by the DFT total energy, 



G«£, 



(16) 



where G and E stand for one of Gxix/Ti an d ^Tix/Ti or 
Gxix/x an d ^Tix/x- This approximation neglects vibra- 
tional contributions (TS ) and the pressure term (pV) 
in G. However, the effect of neglecting these factors is 
expected to be smal L 16 ' 46 



Evaluation of the gas-phase chemical potentials ([15]) 
requires the specification of the deposition temperature 
and the individual partial pressures. The temperature 
and the total pressure inside the CVD chamber can be 
controlled reasonably well during fabrication. Individual 
partial pressures, on the other hand, are not directly ac- 
cessable. However, they can be determined from a rate 
equation that describes the overall CVD process. 

Using the ideal gas approximation (pV — Nk B T), the 
change in individual pressures is described by 



dtPt(t) 



k B T 

VbvD 



Ci Rsuu 



Pi(t) 
Pit) 



R 



cxh 



ViR T 



iX 



, (17) 



where p(t) — JVftW * s the total pressure, Pi(t) denotes 
the partial pressure of species i at time t, and vi is the 
corresponding stoichiometric coefficient (change in num- 
ber of molecules) in the deposition process, Eq. (TTJ). Per 
definition^ stoichiometric coefficients are negative for 
species that are consumed in the reaction and positive 
for reaction products. The value of Ci reflects the con- 
centration of species i in the supply gas. Finally, i? sup 
and i? ox h are the rates at which the gases are supplied 
to and exhausted from the chamber; i?xix is the rate of 
TiX deposition (per stoichiometric formula). 

In steady state, we can eliminate either one of the three 
rates. Keeping the supply and deposition rate as funda- 
mental variables, the steady-state solution {dtp = and 
d t pi = 0) of (flTf is given by 



Pi=P 



Vi ^TiX 



1 + Av rxij 



(18) 



Here, p and pi (without time argument) denote the 
steady-state total (which can be controlled) and partial 
pressures and Av — ^ i>i is the number of molecules that 
are created in each reaction. We have also introduced the 
scaled TiX deposition rate rxix = Rtix/ Rsup- 

Collecting equations pTL2|) . (fTTj) . flT51) . (|T6"j) . and (fTgl) . 
we can express the Gibbs free energy of reaction for de- 
position of excess layers as a function of three variables: 
the temperature T, the pressure p and the scaled TiX 
deposition rate rxix 

G r = G r ({/Xi(T,Pi(r TiX ))}) = G T (T,p,r). (19) 
Here G r denotes one of Gj lX ^ Tl or G^ lX//X . 



III. 5. Limits on the scaled deposition rate 



The scaled deposition rates rxic and rxiN possess crit- 
ical values r T Yc and r"^, i.e., upper values for the pa- 
rameter r that are still compatible with actual growth. 
Deposition of TiX is favorable only if the inequality © 
holds. For rxix > r xiX; the accumulation of reaction 
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products makes the reaction thermodynamically unfa- 
vorable. Thus, in a steady state, r Tix < r^Yx must be 
fulfilled whereas equality corresponds to (dynamic) equi- 
librium. This leads to the conditions 

ffT ic ^MTiCuKYc) + MCH4(4"c) - +4/iHci(4 r ic) ; 

(20a) 

5Ti N ^MTiCk(r^) + 2/x H3 (r^) + ~A* Na 

- 4/* H ci(r?ft), (20b) 

using /itj(r^Yx) = ^(^M^TixD as shorthand. Critical 
values, i.e., upper limits for the scaled deposition rate, 
are identified by solving these equations. 

III. 6. Equilibrium limit 



for the Ti- and X-terminated surfaces, respectively. The 
numbers n\ and n 2 correspond to the numbers of Ti 
layers in the TiX/Ti and the TiX/X slabs, so that 
ni — ri2 = 1. In equilibrium, the difference in free energies 
of reaction is therefore simply the difference in surface 
free energies between the Ti- and X-terminated surface 
At the same time, we emphasize that, 

MTi + Mx — ffTiX ^ MtT + Mx S (26) 

applies for general growth conditions. A nonequilibrium- 
driven process like growth will not generally adjust itself 
to dynamic equilibrium. Out of dynamic equilibrium, a 
prediction of surface terminations in terms of surface free 
energies (employing fix — Mx 8 ) ^ s n0 ^ justified. Instead, 
the free energy of reaction is the proper quantity to use 
to predict surface terminations. 



The AIT-DG method formally resembles the AIT-SE 
when taking the limit in which the gas environment is in 
equilibrium with the surface. Strictly speaking, the AIT- 
DG method expresses the same thermodynamic prefer- 
ence as does a generalization of the AIT-SE method to 
dynamic equilibrium. The dynamic equilibrium is nev- 
ertheless still specified by solving the rate equations for 
the steady-state gas flow. 

We express the differences in reaction free energies for 
excess Ti and excess X surface layers as 

^fTiX/Ti WTiX/X n , aig 

Lr r — (_r r — U-TiX/Ti — A*Ti 

- (G TiX /x - A4 iS ) ■ (21) 

Equilibrium between the gas phase and the surface, 
Eq. (fl"4")) . implies 

= M S rf + M s x urf -M¥f-Mx E ! (22) 

where /iff £ = G T ix/Ti - G T ix and /i^ lf = G Ti x/x - G T ix 
are the chemical potentials of Ti and X at the surface of 
the solid. 

If we also assume equilibrium between the bulk and 
the surface, that is, /i^f* + fi^ rl = /ixi + Mx, and use ([7]), 
we obtain 

Gj 1 - G x = G T ix/Ti — G T ix/x — 5Tix + 2^x S (23) 

as a predictor of preference of equilibrium surface termi- 
nations. 

This equilibrium predictor can be expressed in terms 
of a difference in surface free energies, thus making a 
connection to the AIT-SE^ 

Gj 1 -Gf = 2A ( 7 TiX / Ti - 7 TiX / x ) . (24) 

Here, we use the standard definition of the surface free 
energy 

7 TiX/Ti = ^ (Gxix/Ti - n l9TiX + M x) , (25a) 
7 Tix/x = 1 ( Gtix/x _ „ 25T . X _ Mx ) ; (25b) 



IV. COMPUTATIONAL METHOD 

All DFT calculations are performed with the plane- 
wave code dacapo^ using ultra-soft pseudopotentials^ 
the PW91 exchange-correlation functional; 51 and 400 eV 
plane-wave cutoff. 

The lxl TiX(lll) surfaces are represented by slab ge- 
ometry within a supercell including ~ 10 A of vacuum. 
A 3 x 3 x 1 Monkhorst-Pack k-point sampling^ is used 
and atomic relaxations are performed until interatomic 
forces no longer exceed a value of 0.03 eV/A. The slabs 
contain 12 Ti and 12 X layers for the stoichiometric refer- 
ence surface without excess layer, 13 Ti and 12 X layers 
for the surface with Ti excess layer, and 12 Ti and 13 
X layers for the surface with X excess layer. We have 
checked that the bulk energies per TiC unit, calculated 
as exix = E(N + l) — E(N), are converged to a difference 
of less than 0.03 eV with respect to the slab thickness. 
Here E(N) and E(N + 1) is the energy of a TiX, TiX/Ti, 
or TiX/X slab that contains N and N + 1 Ti layers, that 
is, the bulk energy is independent of the surface termi- 
nation, as it should be. 

For molecules and atoms, we use spin-polarized cal- 
culations when relevant. Total energies are calculated 
within unit cells of size 20 x 20 x 20 A 3 , assuring that 
molecules from different cells do not interact. We use a 
lxlxl Monkhorst-Pack k-point sampling^ and per- 
form atomic relaxations until interatomic forces no longer 
exceed a value of 0.01 eV/A. 



V. RESULTS 

Table |T] lists the calculated atomization energies, bond 
lengths, and bond angles of molecules that are relevant 
for CVD of TiX. For atomization energies, our calculated 
values are in very good agreement for H2 and N2. For 
other molecules, we find larger discrepancies, in partic- 
ular for TiCU. However, since the geometric properties 
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molecule 




[eV/molecule] 


A-Eatom [%] 


b [A] 




A6 [%] 


6 


[°] 




present 


exper.^ 




present 


exper.^ 




present 


exper.^ 


H 2 


4.58 


4.48 


2.2% 


0.75 


0.741 


1.2% 






TiCU 


21.31 


17.84 


19% 


2.18 


2.17 


0.05% 


109.3-7 


109.471 


N 2 


9.63 


9.76 


1.3% 


1.12 


1.098 


2.0% 






CH 4 


18.23 


17.02 


7% 


1.10 


1.087 


1.2% 


109.5 


109.471 


HC1 


4.88 


4.43 


10% 


1.29 


1.275 


1.2% 






a Ref.lH 



TABLE I. Calculated and experimental data for molecules relevant to CVD of TiX. Total energies are required for reference 
for the calculation of chemical potentials. Calculated geometric data, such as bond lengths b and bond angles 8, is in excellent 
agreement with the experimental one. Calculated atomization energies -Batom of N 2 , H 2 , and CH4 differ only slightly from 
experimental values (We use the value of _Eh = —13.6 eV for the H atom in favor of the DFT value to calculate atomization 
energies of molecules that carry H). For Cl-containing molecules, the differences are larger. 




800 900 
T[°C] 

FIG. 4. (Color online) Critical values of the scaled reaction 
rate r^N as functions of temperature at different pressures 
and HC1 concentrations. Red (dark) lines correspond to a 
deposition pressure of pi = 50 mbar, green (light) lines to 
P2 ~ 500 mbar. Solid lines correspond to no HC1 content 
chci.i = in the supply gas, dashed lines to a content of 
1% HC1 chci.i = 0.01. The insert shows the corresponding 
analysis for TiC. 



are in excellent agreement with experiment, we suspect 
that these discrepancies are mainly related to total en- 
ergies of the isolated atoms. We therefore trust total 
energies of the listed molecules as reference energies in 
the calculations of chemical potentials. 

Figure 0] displays the critical scaled deposition rates 
r TiX as functions of temperature for two different pres- 
sures and two different concentrations of HC1 in the sup- 
ply gas. We find that r^"x first increases with increasing 
temperature and then approaches an asymptotic value. 
At fixed temperature, r^jx decreases with increasing to- 



tal pressure. We also find a divergence in log 



crit 
10 'TiX 



(T) 



as the temperature decreases. This divergence can be 
seen most pronounced for TiC in the case where no HC1 
is supplied (divergences of other graphs lie outside the 




bulk TiC 

■ excess Ti - 

■ excess C 



J I L 



J I L 



_L 



log 



-5 

10 r TiC 



FIG. 5. Free energies of reaction G r of CVD bulk TiC (solid), 
excess Ti (dashed), and C excess layers (dashed-dotted) at the 
TiX(lll) surface as a function of the scaled reaction rate rxic- 
We assume a total pressure of p — 50 mbar and a deposition 
temperature of T — 980° C inside the CVD chamber. Red 
(dark) lines correspond to experimental values^ for the sup- 
ply gas compositions, see also main text. In particular, the 
supply gas contains 1% HC1. The green (light) lines corre- 
spond to the case where no HC1 is supplied. 



plotted regime). This divergence arises because lowering 
of the temperature below a certain value causes gJ lX to 
become positive for all choices of rxix- 

TiC(lll). Figure [5] reports the calculated free en- 
ergies of reaction of Ti and C excess layers as well as 
the free energy of reaction of one stoichiometric unit of 
bulk TiC. We plot the free energies of reaction as func- 
tions of the scaled reaction rate rxic at fixed tempera- 
tures and pressures. The scaled reaction rate is limited 
to right by its critical value when dynamic equilibrium is 
reached, see (|20p. We assume experimental values for the 
supply gas composition as stated in Ref. 34. In detail, 
we have cticm = 0.04, cch 4 = 0.07, chci = 0.01, and 
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FIG. 6. Free energies of reaction G r of CVD bulk TiN, excess 
Ti and excess N layers at the TiX(lll) surface as a function of 
the scaled reaction rate r-riN- In the top panel, experimental 
values^ for the supply gas compositions from are used, see 
also main text. Temperature and pressure are chosen as in 
figure [5] The termination critically depends on the the scaled 
deposition rate. The bottom panel illustrates the effect of 
varying deposition parameters. For bulk, excess Ti, and ex- 
cess N, the same line style as in the upper panel applies. Blue 
(dark) lines correspond to an HC1 concentration of 2% in the 
supply gas. Green (light) lines correspond to the case where 
the N2 concentration is decreased to 1% and the temperature 
is raised to T = 1200° C. In both cases, the changes are as- 
sumed to be balanced by the H2 concentration. An increased 
HC1 concentration leads to a N-terminated surface, indepen- 
dently the value of rxiN- A decreased N2 concentration favors 
Ti-termination. 



CH 2 = 1 - Y,i X i- 

We find that Gj lC ^ C > Gj lC ^ Tl in the entire range of 
the scaled deposition rate rxic- Furthermore, the free en- 
ergy of reaction associated with a C-terminated TiC(lll) 
surface is positive. Thus, under the considered exper- 
imental circumstances, the TiC(lll) surface will be Ti 
terminated. 

We have also tested the consequences of varying depo- 
sition parameters. With the chosen set of supply gases, 



it was not possible to identify a set of deposition param- 
eters that could lead to a favored C-termination of the 
TiC(lll) surface. 

TiN (111). Figure [6] reports the calculated free energies 
of reaction of Ti and N excess layers as well as the free 
energy of reaction of one stoichiometric unit of bulk TiN. 
In the upper panel, we assume experimental values for 
the supply gas composition as stated in Ref. 1371 In detail, 
we have xticw = 0.09, xn 2 = 0.5 and xh 2 = 1 — Xi. 

We find that the free energies of reactions of both ex- 
cess Ti and excess N can be negative simultaneously. Far 
from dynamic equilibrium, that is, for rxiN << ^t?N' we 
have Gj lN ^ N > Gj lN ^ Tl . In this regime, a Ti-terminated 
surface is more favorable. Close to dynamic equilibrium, 
that is, for tun ~ '"tin ' we have Gj lN ^ N < G™^ Tl . 
Furthermore, £^F lN / Tl is positive there. In this regime, a 
N-terminated surface is more favorable. 

In the lower panel we illustrate the effects of vary- 
ing the deposition parameters. Assuming an increased 
HC1 content in the supply gas, balanced by the H2 con- 
tent, we find that Gj lN ^ Tl turns positive over the whole 
range of the scaled reaction rate. Thus, we predict a N- 
terminated TiN(lll) surface. Conversely, a decrease of 
the N2 concentration in the supply gas (balanced by H2) 
and a simultaneous increase of the deposition tempera- 
ture, to T = 1200°C results into the favorization of a 
Ti-terminated surface. 

In summary, we find that, under the conditions speci- 
fied in Ref. |34], the as-deposited TiC(lll) surface will be 
Ti terminated. Forcing a C-tcrminated surface by varia- 
tion of process parameters within the specified supply gas 
can not be achieved. Under the conditions specified in 
Ref. [37l the TiN(lll) surface can, in principle, be either 
Ti- or N-terminated, depending on the ratio ttin between 
the deposition and supply rate. By introducing a slight 
amount of HC1 into the supply g N-terminated sur- 
face will be forced. With another choice of deposition 
parameters a Ti-terminated surface will evolve. 



VI. DISCUSSION 

The AIT-DG method that we have proposed in Sec. lIIII 
and illustrated in Sec. [V] is designed to describe surface 
terminations as a function of the growth environment. 
Our predictions are based on calculations of the free en- 
ergies of reaction G r . We motivate this criterion by the 
physical principle that all systems strive for the states 
of low energy, given relevant boundary conditions. The 
method applies in dynamic equilibrium but also in sce- 
narios where (dynamic) equilibrium is not maintained 
but the system still settles into a steady state. The 
method assumes that kinetic effects do not lock the sys- 
tem into thermodynamically unstable morphologies. It 
assumes in essence that a high deposition temperature 
causes a continuous annealing of the emerging structure. 
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VI. 1. Model limitations 

Kinetic-barrier effects certainly exist and may prevent 
some systems from actually reaching the structure and 
morphology with the lowest thermodynamic energy^ 
Similar to the AIT-SE method, our method neglects 
these barrier effects. In principle, this approximation lim- 
its the applicability of our method to high-temperature 
growth where barriers become less important. Never- 
theless, kinetic effects can not be expected to alter the 
predictions in cases similar to TiC. This is because one of 
the free energies of reaction is negative (Ti excess) while 
the other is positive (C excess). 

Another obvious limitation is that of the ideal-gas ap- 
proximation from which we infer the temperature and 
pressure dependence of the molecular chemical poten- 
tials. Formally this approximation limits our description 
to low-pressure scenarios, but the approach could easily 
be generalized by a more refined description of the gas 
thermodynamics. 

VI. 2. Growth and rate equations 

Our method is capable of describing surface termina- 
tions as they evolve during a growth process. This is fun- 
damentally different from predicting the termination of a 
surface that is in equilibrium with the environment. As 
illustrated for TiN in the top panel of figure [H the emerg- 
ing surface termination depends on the the scaled reac- 
tion rate, r-rix = Rnx/ Rsup- In such cases, our method 
should have clear advantages over the AIT-SE method. 
Use of our method requires only that measurements of 
the scaled reaction rate are made available^ 

We are dealing with an open system that has both a 
gas flow and a consumption of gases during deposition. 
In a closed system, where the number of gas-phase species 
is fixed, the surface-gas system will approach equilibrium 
(unless there are important kinetic barriers). 28 Similarily, 
the growing surface and the gas environment could be in 
a state of dynamic equilibrium. In this case, the analysis 
in Sec. IIII.6I shows how the AIT-SE formalism could be 
extended to cover a description of the deposition. This 
can even work for materials that do not contain con- 
stituents A for which there exist molecular counterparts 
A n 4£ Using the free surface energy as equilibrium pre- 
dictor, the chemical potentials of atoms in the gas can be 
associated as an appropriate difference between molecu- 
lar chemical potentials. An explicit example for TiC in 
an CH4-TiCi4-H2 environment is given in Eq. (fUZj) . How- 
ever, since the chemical potentials of atoms now depend 
on partial pressures of several gas-phase species, there is 
still a need to determine all relevant partial pressures. 
The rate-equation approach presented here supplies the 
information that is necessary for application of either the 
AIT-SE or the more general AIT-DG description of de- 
position growth^ 

We emphasize that growing systems, in general, do not 



need to reach dynamic equilibrium, but may be in any 
other steady state (if steady-state is reached at all). If 
dynamic equilibrium was always maintained, the scaled 
reaction rate rxix would always assume its critical value, 
that is, a maximum. The absolute deposition rate could 
then be increased to arbitrarily large values, simply by 
increasing the supply rate at a fixed total pressure. We 
believe that the assumption of dynamic equilibrium is 
too optimistic for CVD systems, and we would not gen- 
erally trust dynamic-equilibrium descriptions. The here- 
proposed theory of deposition growth (AIT-DG) applies 
also when the deposition and supply rates are indepen- 
dent variables. 



VI. 3. Innovation potential 

Using the free energy of reaction for the formation of 
excess layers on a specified surface we have shown for 
TiX(lll) how surface terminations can be understood 
from knowledge of the details of the CVD gas supply. 
For TiN, we have also illustrated that the method can be 
used to guide the growth towards a desired surface com- 
position. That is, we have shown that we can modify 
the deposition environment (supply gas composition, de- 
position temperature, and deposition pressure) to allow 
either of the two surface terminations to emerge. 

From this point of view our suggested method has a 
potential for accelerating innovation^ In the case of 
TiX(lll), the potential to design the surface termina- 
tion may seem trivial. The change from a Ti- to a N- 
terminated TiN(lll) surface that follows from raising the 
HC1 concentration in the supply gas could likely also be 
obtained by, terminating the deposition process with an 
replacement of the TiCl4-N2-H2 supply gas with a pure 
N2 supply gas. We argue, however, that our method is 
valuable for characterizing more complex materials with 
a larger variety of possible surface terminations. This 
is particularly true when we extend the application to 
understand binding at interfaces that arise in deposition 
growth^ for example, in the case of thin-film alumina 
on TiC^ 



VII. SUMMARY AND CONCLUSIONS 

We present a method based on ab-initio calculations 
supported by thermochemical data to predict termina- 
tions of surfaces as they emerge during high-temperature 
deposition growth from a multicomponent gas, for exam- 
ple chemical vapor deposition. The method describes sce- 
narios that are located between a static equilibrium and 
the fully non-equilibrium regime that requires a complete 
kinetic description. It relies on the basic thermodynami- 
cal principle that all systems strive for low-energy states 
subject to statistics and other boundary conditions. Rate 
equations and associated steady-state solutions play a 
central role in the calculation of free energies of reaction. 
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We use the latter to determine the chemical composition 
of the outermost surface layer of the growing surface. 

We illustrate the approach for TiX(lll) (X = C or 
N). For TiC, we predict a Ti-terminated surface when 
fabricated under the conditions stated in Ref. HU- For 
TiN, our predictions based on experimental conditions 
reported in Ref. H3 are not so clear. We find that the 
termination depends on the ratio between the reaction 
rate and the rate at which the gas is supplied to the re- 
action chamber. 57 In or close to dynamic equilibrium, we 
predict a N-terminated TiN surface. Departing from dy- 
namic equilibrium will result in a Ti-terminated surface. 
We also suggest deposition parameters for which Ti- or 
N-terminated surfaces can be achieved independent of 
reaction and supply rates. 

We also compare our method with the AIT-SE 
method. 15 ' 16 We show that the predictions of our method 
agree with those of the AIT-SE method in the limit 
where (dynamic) equilibrium is maintained, but also 
point out that this condition is not necessarily applicable 
to growth. In a closed system (absence a gas flow), there 
is an affinity to reach equilibrium. In an open system 
(with the gas flow turned on), one could expect a cor- 
responding affinity to reach dynamic equilibrium. How- 
ever, as discussed in Sec. EH this expectation is likely too 
optimistic for general CVD growth conditions. In such 
cases, our method has clear advantages since it does not 
a priori implement equilibrium constraints. 

The choice of materials, in particular TiC, also ex- 
emplifies the broader applicability of our method as 
compared to a strict implementation of the AIT-SE 
method^ Our method is not formally limited to pre- 
dictions of oxide surface stability (or as in extensions of 
AIT-SE, to A TO B„ compound surfaces for which there is 
an Afc or Bfc counterpart in the gas phase). 

Finally the capability of predicting different surface 
terminations as a function of the deposition conditions 
illustrate the predictive power of our method. It moti- 
vates continued development of ab-initio thermodynam- 
ics, also in a role to guide experimental optimization of 
surface and interfacial structures. 
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Appendix A: Gibbs free energy of formation as 
predictor of surface terminations. 

Chemical reaction theory provides the basis to relate 
the probabilities Pa and Pb for A and B surface termina- 
tion to free energies of reaction G r . The surface transfor- 
mations described by (J5J) can be cast into the following 
rate equations for the probability of growing either an A- 
or B-terminated surface, 

d t P A = -(Tl + T 1 / )P A + {Tj + Ti 1 )P B (Ala) 

d t PB = (ri + rf )p a - (rj + r^)p B . (Aib) 

The steady-state solution is P A /P B = (rj- + r£ 7 )/(r£ + 
Ty 7 ) where rj^p are forward (f) and backward (b) reac- 
tion rates identified in reaction (|5aj) and (I5bp . These rates 
can be expressed as the products of the corresponding mi- 
croscopic rate constants ^Y/iP 1 anc ^ °^ an environment- 

specific product Q {fM = II i [Xi] !/ "' b " *, Ref . [H and [U. 
Here pQ] denotes the concentration of molecule Xj and 

v Yfb\ \ 1S ^ nc (positively counted) stoichiometric coeffi- 
cient as reactant (f) and product (b) in reaction (|5a|) or 
([5b|l. We note that Q{f,b} is strictly speaking a product 
of activities from all parts in the reactions (|5a|) and (|5b[) , 
but for surfaces (or solid) the activity can be approxi- 
mated by unity. 

The general relation (3G r = — InK + lnQ, with j3 
the inverse temperature, relates G r to the equilibrium 
constant K — kr/kf, and Q — Qb/Qf- We have 

Q {/b} } = T \Yb} l k ifM and find exp [G? - Gf]) = 

Dynamic equilibrium is characterized by Gf + Gj 3 = 
or equivalently T^ry / (TItI 1 ) = 1, from which it follows 
that 

P aAPb | dyn . oq . = exp (-/3 [Gf - Gf] /2) . (A2) 

Away from dynamic equilibrium, this relation may 
not hold exactly. However, whithin a broad range 
of noncquilibrium conditions we may still approx- 
imate the arithmetric mean [{x, u)am = (x + 
y)/2] in Pa/Pb = {Tj, T^am/^, rf) AM by the 
geometric mean [(x,i/)gm = (xy) 1 ^ 2 ] expressing 
exp (-/3 [Gf — Gf] /2) = <r}, T^gm/^, T y ) GM . 
This suggests that the evaluation ([6]) remains a good ap- 
proximate measure of relative surface termination also 
out of equilibrium. 
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